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SYNOPSIS 

Digital Sim-ulation of Beas Satlej System ” - 
a thesis sutmitted in partial ftilfilment of the requirements 
for the Degree of Master of Technology hy Vljendra Singh to 
the Department of Civil Engineering, Indian Institute of 
Technology, Kanpur, November 19'^8* 

Water resources systems are .generally large and 
complex. They consist of multiple units and serve multiple 
purposes- The multiple purposes are not wholly complementary. 
Several general or problem specific simulation models and com- 
puter programmes have been developed for analysing water re- 
sources systems. SIMYLD II is a computer programme developed 
by Texas Water Development Board for simulating the hyd-rologic 
operation of a system of interconnected reservoirs within 
a basin or a multibasin water resources system. The study 
consists in implementing the SIMYLD II programme^ validating 
it with available data and adapting it for the operation of 
Bhakra Beas system. 

The original SIMYLD II programme was implemented 
in IBM 7044-1401 system at I.I.T. Kanpur and was adapted to 
meet the requirements of Bhakra Beas system. 

Using 13 years of historical data, the Bhakra 
Beas system was simulated using the modified SIMYLD II model. 
The criteria for defining wet, average and dry years; the 
rule curves for operation in i^ret, average and dry years; and 
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the relative weight ages for meeting different demands and 
for maintaining the rule curves are derived from simulation 
analysis. Results indicate that by using these criteria, 
the benefit from the operation of the system can be greatly 
increased. Further improvement of the model is also possible 



1. IHTROnJCKON 


1.1 General 

1.1.1 Water resources planning 

Water resources systems are generally large and complex. 
IThey consists of multiple units ard they serve multiple pur- 
poses. The multiple purposes are not wholly complementary 
Water resources systems are designed to serve several socio- 
economic ohjectives like national and regional economic deve- 
lopment, income distribution, preservation and athancement of Envi 
nonmental quality and social well being (Maass et al., 1962, 
United States Water Resources Council, 1973). They are also 
affected by economic and hydrologic uncertainties. 

The conventional practice in planning, design and 
operation of water resources projects has been to consider a 
set of demands and to satisfy them by river basin development. 
Such a development includes the purposes it is to serve, the 
physical means for meeting these purposes, the sizes of needed 
facilities and the levels of output-, Based on experience, a 
few alternatives are proposed, analysed and evaluated through 
incremental analysis before a final design is adopted. 

Planning can be defined as the orderly consideration of 
a project from the original statement df purpose through the 
evaluation of alternatives to the final decision on a course 
of action. It is the basis for decision to proceed with a 
proposed project and is clearly the most important aspect of 
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the total engineering for the project. The planning for an 
entire rirer hasin involves a much more complex planning 
effort than that for a single project, hut the diffioulties 
in arriving at the correct decision may he jnst as great for 
the smaller project. 

Huf Schmidt considers the field level planning of water 
resources systems to consist of four hierarchical steps, viz. 
(Gold, 1975)- 

i) Definition of project objectives and agreement of these 
with policy makers,* 

ii) Securing staff and scheduling their tasks J 

ill) Data assembly on the physical variables (hydrology, 
water quality), human requirements (water supply, land use, 
recreation),' and economic quantities (construction and opera- 
tion costs, explicit benefits),* and 

iv) formulation of system design via (a) screening of alter- 
native configurations and (b) detailed analysis of remaining 
components. 

Planning may deal with modification to and operation of 
an existing system, or the design of an essentially new system 
involving the operating variables also. The latter is naturally 
more complicated. The two parts of step iv imply that 
initially a large number of alternatives are to be considered, 
and so relatively crude models are used to essentially elimi- 
nates trivial and inferior alternatives and isolate better 
alternatives. This is generally done by programming models. 
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Onee better alternatives are identified, it may be possible 
to consider more detailed mcfdels of the system and arrive at 
the ” best" one. This may involve stochastic programming and 
simnlation models. 

Wiener (1972) considers the role of planner and system 
analyst in the socio-economic development process in an xinder- 
developed conntiy where water is used as a critical input for 
economic development and progress. He deems it necessary to 
consider the implementation phase also in planning particularly 
in the case of underdeveloped countries. He also indicates the 
limitations and opportunities in the application of system 
analysis and operations research techniq_ues in water resources 
planning. 

1.1.2 Simulation analysis 

simulate means to duplicate the essence of a system 
or activity without attaining reality itself. Simulation has 
been used traditionally in engineering. The use of conceptual 
system models, scale models, analogues and laboratory experi- 
mentation are but some of the general simulation techniques 
traditionally used in engineering. Simulation has been used 
for a number of purposes including the analysis of the system 
to estimate the parameters and the behaviour of the system that 
is existing or is yet to be* the effect of the environment on 
the system design^ the demonstration of the performance of a 
new complex system, and for giving training in the control of 
complex systems. 
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Digital simulation is the numerical simulation of 
the process in a digital computer. A behaviour model of 
the physical components of a water resources system is for- 
mulated in terms of the components ^parame tens, variables 
and relationships among them to study the processes as they 
■'evolre in time through the several components of the system when 
subjected to a series of hydrologic inputs and human 
interventions in terms of design, construction and operation 
of the system. When the mathematical model for a process has 
been decided upon, the various elements of the process can 
be represented on the computer so that outputs from one part 
of a system constitute inputs to one or more elements of 
the system or to itself. The system is simulated for a set 
of design and operating parameters, and the effect of changes 
in these parameters on the system response is Investigated. 

The simulation of any system involve generally a 
number of steps which can be represented as follows(Fig.l.l) ? 


i) formulation of the problem The problem is to be 
formulated in analytic terms. This involves the definition 
of the objective, their priorities, the approach to optimiza- 
tion, the identification of the system and its environment 
and the structure of the system. 

ii) Analysis of data ^ The historical data for the inputs, 
the components of the system, and the outputs are collected. 
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ill) Formulation of system model Initially component 
models are formulated; from the data of the previous step, 
parameter values are estimated for the component models; and 
the models are validated- A flow diagram for the system 
indicating the sequence of the modifications of the inputs 
hy the system is prepared. 

iv) Estimation of parameters and validation of system 
model : The parameters of the system are estimated, the 
system is analysed or numerically simulated and from the 
responses of the system for hypothetical or historical inputs 
as per record and from the results of simulation, the formu- 
lated system model is validated. 

v) Sampling for design variables *. A set of feasible and 
preferred design and operating parameters are chosen. 

Vi) Simulation of the process '. The system is digitally 
simulated using historical or generated data as the case may 
be . 

vii) System utility The measure of system response in 
achieving the goals for the system parameters assumed in each 
case is evaluated- 

The results of simulation are also analysed to deter- 
mine whether the system design can be improved with reference 
to the achievement of the goals. The parameter values are 
modified, and steps vi and vii are repeated until the "best” 
of the alternatives considered is identified- 
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Yiii) ReiDO?'ting of J'e suits ffhe results of the study, the 
conclusions drawn and the recommended system design are 
reported in the final stop for implementation. 

1.1-3 Planning models of IlOB 

SeTeral general or problem specific simulation models 
and computer programmes have been developed for analysing 
water resources systems. Important references include Maass 
et. al. (1962), Hufschmidt and Piering (1966), lucia et al 
( 1971 ), U.S. Army Corps of Engineers Hydrologic Engineering 
Centre (1968), and Texas ¥ater development Board (1972, 1972a, 
t974-). The general approach to simulation may be based on 
heuristic criteria as in the case of Maass et.al.(l962), 
Hufschmidt and Piering (1966) or may use an imbedded optimi- 
sation technique as in the case of Texas ¥ater Resources 
Board (1972), 

Texas ¥ater development Board (T¥BB) "fesgsn. in' 1967 .■ 
long range programme of applied research in water resources 
system simulation and optimization. The objective was to 
develop a set of generalised computer oriented planning tools 
for use in detailed planning, design and management of water 
resources systems such as the Texas Water Systen as proposed 
in the Texas Water Plan. Project II refers to the second 
phase of a three phase research project leading towards the 
development of a computer oriented planning system for use 
in planning of large, multibasin systems of reservoirs and 
connecting river reaches and pumped canals. The six planning 
stejjs of project II are shown in Pig. 1.2. They include ^ 



USED 



JNDEMTtnCATION OF 

OBJECTIVEi, AND GOAL5 
AND 

DEVELOPMENT OF 
STOCHASTIC DATA BASE 


STEPS ONE 
AND TWC TEN ^ 
DEVELOP Pk-- .:.P,hMMES 
OBJECTIVES, ANC 
DATA ■ BASE 


INDIVIDUAL BASIN AND 
MULTI BASIN 
FIRMYIELD 
ANALYSIS 

RESERVOIR sizing;) 


SINVLO-H 


multibasin 

PLAN improvement 
AT 

PRESPECIFIED 

POINTS 

ON 

DEMAND -BUILDUP CURVE 


fSTEP FOUR 
■A PLAN 
I IMPROVE WENT 


MULT! BASIN 
.-■-ilNG PLAN 

OPTIMIZATION 
TO FIND 

w- E y p ECTEDi- COST 
STAGING PLAN 


jRlC'AL AND 

HASTIC 

BASE 


j 'SENSITIVITY ANALYSIS 

. • INIMUM-EXPECTED-C OST 

I PLAN 

, IaND viable ALTERNATIVES 

! ; 

1 t SELECT PINAL 
I implementation plan 


STEP SIX SIM III 
FINAL PLAN AND 
SELECTION AL-II 


1-2 f KOCEDURAL PLANNING STEPS FOR PROJECTU 





9 


Step 1 I Identification of objectives and goals - This step 
outlines the problem to be solved, specifies in general the 
magnitude and location of demands to be met and the priori- 
ties associated with meeting each of the specified demands, 
identifies the sources of the water to be considered, and 
identifies the criterion to be used in optimising the 
selection of an implementation plan- 

Step 2 Data base development - In this step both histori- 
cal and stochastically generated hydrologic data sets are 
developed concurrent with the economic and physical data- 

Step 3 • Plan development - This step analyses each of the 
individual river basins to determine their firm yield 
characteristics and the location of possible attractive basin 
import and export points. 

Step 4 • Plan improvement - This step improves those plans 
such that minimum - cost plans are developed at each of 
several points on a prespecified demand - build up curve. 

Step 5 Plan optimisation - It optimises the staging plan 
over the time period when demands for water are increasing. 
This implies finding the minimum - expected cost plan that 
meets the demands specified with an optimal level of shortage 

Step 6 I Pinal plan selection -■ This is the final step in 
the planning procedure and involves testing the sensitivity 
of the cost and physical response of the simulated prototype 
to variation in all of the important parameters used by the 
models. 
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1.1.4 SIMYID II progiraiiiirie 

SIMYID II Is a comp-ater programnie designed to simulate 
the hydrologic operation of a system of interconnected reser- 
voirs within a hasin or a multibasin water resources system. 
The inputs to the model are physical description of the system 
to simulate operating criteria and monthly inflows, monthly 
demands and monthly evaporation rates. The output consists 
of monthly end of month reservoir storages, monthly flows in 
the system^ s river reaches and canals, and annual end period 
of simulation summaries. The main attributes of the model are 
simplicity in set up and use, generality in application, and 
speed of computation. Other models are more comprehensive 
and naturally much more complicated. Since SIMYII) II is 
used in earlier steps in decision making than the other models 
(lig. 1.2) and it is also useful in identifying surplus and 
deficits, it is proposed to implement SIMYII) II at IBM 7044- 
1401 system at 1. 1. T. Kanpur and test it for field data of a 
river basin. 

1*2 Statement oflSroblem 

It is proposed to implement SIMYID II programme deve- 
loped by the Ti/ICB, validate the programme using the test data 
given in the TWDB report publication on the model ( TT/i/DB, 1972 )j 
adapt the simulation programme if necessary to meet the 
requirement of water resourdes systems in India jand test and 
validate the adapted model using some data for an Indian 


river basin. 



11 


1.3 Objective of the Study 

. The Tuajor objective of the study is to implement and 
adapt the SIMYID II programme of nCB to suit requirement of 
a river basin in India, identify the limitations and capa.bi- 
lities of such a model and hence develop knowhow for more 
comprehensive development planning models. 

1.4 Scope of the Study 

The scope of the study is limited by i) available data 
for real systems’ ii) time availabilityfof the study," iii) 
limited interactions with fi eld engineers and planners^ and 
iv) limitations in the IBM 7044-14Q1 computer system at I.I.T. 
Kanpur • 

Hence the scope is limited to the following 

i) Because of the limited capacity of the IBM 7044-1401 
system, simulation is limited to 15 nodes, 20 links and 13 
ye af s . 

ii ) The data were available for Bhakra-Beas system for a 
number of years and because of the limitation in time, com- 
puter capacity and computer availability, the studies are 
limited to the following . 

a) While Bhakra Beas Management Board (BBMB) adopt a ten 
daily operation of the system, only monthly operation is 
considered in this study' 

b ) Only 13 years of data values are used in this study; 

c) Because of limitation of the data availability and 
non avoidable lumping of demands and supplies in the model. 
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certain asstimptions were made with reference to allocation of 
demands for different sources and for the lumping of demands 
at nodes,* 

d) The demands are defined "by BBMB for dry, dependable and 
average years. The inflows in the system may be dry, arerage 
or wet. Certain reasonable assumptions were made to link the 
demands to the inflows and these are considered in great detail 
in Chapter and 

e) Because of the limitations of the study^ model needs 
further modification as indicated in Chapter 5* 

1-5 Significance of the Study 

As the remaining available uncommitted supplies of water 
and land resources diminish and demands for them increase, the 
objectives of water resources planning broaden "the physical 
facilities required become more complex, and the limitations 
under which they must be implemented become more stringent. 
There exists an urgent need to develop techniques which can 
enhance the capability of the planners to make- an intelligent 
and comprehensive evaluation of alternatives. Because costs of 
construction, operation, and maintenance of water resources 
facilities are likely to be large, a means, must be found for 
analysing alternative solutions to water problem. Planners 
have been turning to sophisticated mathematical techniques 
applied on digital computers of increasing speed and accuracy. 
This study will lead to a better understanding of the problems 
in managing our scarce national water resource and towards a 
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iDetter tool in the design of optimum water resources systems 
in India. 

1.6 Organization of the Study 

The study is reported in the following sequence *. 

i) Chapter 2 discusses the River hasin simulation Model 
SIMYIh II and indicates the capabilities and limitations of 
the original programme) 

ii) Implementation of programme SlitYLI) II is discussed in 
Chapter 3- Modifications, additions and alterations in the 
programme on IBM 7044-1401 systems to suit Indian conditions are 
also described,' 

iii) Chapter 4 discusses the use of the adaptations of 
SIMYII) II in the simulation of Beas-Sutlej System. System 
description and input data required by the model are describ- 
ed briefly. A comparison is made in study of planned operation 
and improvements in operation. Results obtained and tentative 
conclusions are also discussed) and 

iv) Summary, conclusions and suggestions for future study 
are presented in Chapter 5. 
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2. EIVER BASIH SIMUMTIOF M03DEL, SIMYIB II 

2.1 Introduction 

SIMYLD II is a computer programme dereloped tj TTOB 
designed to simulate the operation of the system of inter- 
connected reservoirs within a hasin or a multihasin water 
resources system. The description of the model and its various 
details have been taken from several publications of the TTOB 
(particularly TTOB, 1972). 

2.2 Model Description 
2.2.1 Purpose 

The purpose of SIMYLD-II is to provide the water 
resource planner with a tool for analysing water storage and 
water transfer within a multireservoir or multibasin system. 

The model has the following uses to the planner 

i) It simulates the movement of water in a system of 
reservoirs, rivers, and conduits on a monthly basis while 
trying to meet a set of specified demands in a given order 

of priority. If a shortage occurs (i.e., not all demands 

can be met for a particular time period) during the operation, 
they are spatially located at the lowest priority demand 
nodes' 

ii) It determines the firm yield of a reservoir within a 
water resources system. Pirm yield is defined as the maximum 
demand at a reservoir that can be met with ’acceptable’ 
shortages^ 



15 


iii) The model is designed to proTide the user with flexibi- 
lity in selecting operating rules for each reservoir. The 
operating rules are formulated as the percentage of the 
reservoir capacity (either total or conservation) that. is 
desired to he held in storage at the end of each month. In 
addition to it, priority ranking, used to determine the allo- 
cation of water between meeting demands and maintaining storage 
is assigned to each storage reservoirj 

iv) The model can analyse either static or d3rnamiic system 
operation, in that both constant or time-variable demands can 
be analysed^ and 

v) The user can analyse the operation of the system under 
the expected ultimate demands for any selected hydrologic 
sequence . 

The model is designed to simulate both small scale 
systems, such as two or three reservoirs within one river 
basin, and large scale systems, such as the proposed Texas 
Water System. 

2.2.2 Concepts 

The concept behind SIMYLD-II is that the physical water 
resource system can be transformed into a capacitated network 
flow problem. While making this transformation, the physical 
elements of the real system are represented as a combination 
of two possible network components — nodes and links. 
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As the nomenclat-ure implies, a node is a connection 
and/or "branching point within the network. Therefore, a node 
is similar to a reservoir or non-storage junction (i.e. canal 
junctions, major river intersections, etc.) in the physical 
system. Along with this, a node is a netv/ork component which 
is considered to have the capacity to store a finite and 
bounded amount of the material moving in the network. In the 
case of SIMYIiD-II, reseivoirs are represented by nodes which 
have a storage capacity as well as the ability to serve as 
branching points. A nonstorage capacitated junction is treat- 
ed similar to a capacitated junction (reservoirs) except that 
its storage capacity is always zero. Demands placed on the 
system must be located at nodal points. 

The transfer of water among the various network nodes 
is fulfilled by transfer components called links. A link is a 
river reach, canal, or closed conduit with a specified direc- 
tion of flow and a fixed maximum and minimum capacity. The 
physical system and its basic time step operation, in this 
case a month, is formulated as the network flow problem. The 
network flow problem is nothing more than a mathematical 
representation of the physical system. This mathematical 
representation is accomplished as follows ’. 

i) Reservoirs and nonstorage junctions are represented 
by nodes’ 

ii) River reaches and pumped canals are represented by 
links,' and 
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iii) Additional information needed to descrilDe the system, 
such as inflows, demands, spill points, and starting conditions 
are specified hy the user. 

An initial step in the application of SIMTLD-II is 
the construction of the node-link diagram descnihing the 
physical system. In designing the node-link diagram, the 
physical system elements are represented hy diagram elements 
in the following manner ' 

(i') Reservoirs are represented hy triangles j , 

(ii) Non-storage junctions or branching points are represented 
by circles 

(iii) River reaches are represented by dashe -d lines showing 
the directions of flow ,’ and 

(iv) Canals or closed conduits are represented hy se d - i - d - 
lines showing the direction of flow* 

In order to make the node-link diagram conform to the 
requirements of SIMYLD-II, the following rules for designing 
the diagram should be followed \ 

i) Water can enter or leave the system only at node points 
(either storage or non-storage nodes). Inflows C over land 
flow, tributaries, etc-), link losses ( evaporation, seepage, 
etc.) and demands can be lumped at the closest node. If more 
detail is required, additional non-storage nodes can be 
inserted at critical locations ,* 
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ii) The numbering system used to describe the nodes con- 
sists of numbering all nodal reservoirs consecutively 
followed by numbering all non-storage nodes ‘ and 

iii) Number all river reaches consecutively followed by 
all canals and/or conduits in the same manner as described 
above. Figure 2.1 shows this node arc configuration formu- 
lated as a network flow problem. 

Once the node-link diagram is complete* all nodes, 
links, and other information must be described in terms 
of directed, capacitated arcs and nodes. Within the typical 
network there are seven types of nodes. These are I 

i) Reservoir nodes * 

ii) Non-storage nodes 

iii) Initial storage and inflow nodes * 

iv) Demand nodes j 

v) Spill nodes ,* 

vi ) Final storage nodes * ahd 

vii) Net balance nodes 

Connecting these seven types of nodes are seven types of 
arcs- Flows in these arcs are constrained to be within 
specified upper and lower limits. The seven types of arcs 
are I 

i) Physical system link arcs ( river reaches, pumped 
canals, etc.) ' 

ii) Initial storage and inflow arcs 
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iii) End of montli desired storage (operating rules)* 

iv) Balance of final storage and maximum reservoir 
capacity arcs 

v) Demand arcs J 

vi) Spill arcs j and 

vii) Net "balance arcs. 

Table 2.1. shows the arc types and their upper and 
lower bound constraints- The total number of arcs in any 
network is given by * 

ARCS = + 3(%) +Ng +4 (2.l) 

vdiere 

Nt = number of river reaches and canals, 

-Lj , ■ 

= number of reservoirs and non-storage junctions, 
= number of demand nodes, 

Ng = number of spill nodes, and 

4 = number of balance arcs. 

Cost per unit of flow is associated with each arc 
in the network. These unit cost coefficients are used to 
find the minimum cost solution to the network flow problem. 
As input to the model, the user selects priorities for meet- 
ing demands and satisfying final end-of month storage 
requirements in the reservoirs. These priorities are then 
converted into the above mentioned costs by the programme. 
Demand arc and desired storage arc costs are expressed as 
negative costs ( analogous to benefits). The more negative 
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[TABLE 2.1 ARC [TYPES MD BEEINITIOyS OF THEIR UPPER ANB 
LOMR BOURRS CT!/Jl)B, 1972 ) ’ 



Arc Type 

Lower Bound 

Upper Bound 

1. 

Physical system link 



B. • 

River Reach 

Minimum River Capa- 
city+ 

(User Specified) 

Mazimum River 
Capacity 

b. 

Canal, Pipe line 

Minimum Canal Capa- 
city +• 

(User Specified) 

Mazimum Canal 
Capacity 

2. 

Initial Storage 
and Inflow 

Previous End-of- 
Month Storage Plus 
Current Monthly 
Inflows 

Previous End-of- 
Month Storage 

Plus Current 
Monthly Inflows 

3. 

Einal Desired 

Storage 

Reservoir Minimum 
Pool 

(User Specified) 

Percent of mazimum; 
Capacity Desired : 
(monthly operat- 
ing rules) 

4- 

Final Storage 

Balance 

Zero 

Balance Between ! 

Mazimum Capacity 
and Upper Bound 
of 3 . 

5 . 

Demand Arc 

Zero 

Demand at Node. 

6 . 

Spill Arc 

3e ro 

System Capacity 
Multiplied by 10 i 

7. 

Ret Balance 



B. • 

Total Initial Stora- 
ge Plus Inflows 

I Initial Storages 

1 Initial Storages 

b. 

Total Final 

Storage Arc, 

I Final Storage Balan-X Final Storage 
ce Plus X Final De- Balance plus 
sired Storage £ Final Desired 

Storage 

c • 

Total Demands 

Zero 

£ Demands 

d. 

Total Spills 

Zero 

£ Spills 


Hh 

ShoiAld Be zero 01111633 mininmin flow is required. However a 
minim-um flow requirement may cause infeasible solutions. 
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the riumher, the higher the priority for meeting the upper ; 
constraint ( demand or & borage). Physical system link and 
spill costs, on the other hand, are positive costs. 'The 
effect is to meet demands and desired storage in the order 
of the priorities while minimising the canal pumping and 
spills from the system. It may he noted that high rank in-l 
dicat es ih:.: lower priority and vice-versa. 

2.2.3 Assumptions 

i) Evaporation losses for all reservoirs are calculated 
by the product of the monthly evaporation rate and the aver 
monthly reservoir surface area ,* 

ii) Demands for water are known for the month being 
simulated ' 

iii) Unregulated inflows to the system are known for' the : 
month being simulated ,* 

iv) Reservoir storage contents are allowed to fluctuate : 
between the maximum and minimum capacities specified by 
the user j 

v) Spills occur only at specified nodes and are the mos; 

expensive alternative ) and : 

vi) The flow in all links ranges between the maximum and; 
minimum capacity specified by the user* 

2.2.4 Steps 

The procedure adopted in the programme makes use of i 
the following four steps in moving from a known set of stat€ 
variables at the beginning of a time step to the solution ; 
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for the required set of state 'vax'iables at the end of the 
time step. The four solution steps are summarized as follo^vs.’ 

i) The present status of the network is evaluated and 
all system elements are given an appropriate parametric 
description , 

ii) All specified h^rdraulic and hydrologic inputs and 
demands are accounted for, and the mass halance for the 
entire network system is determined. Bounds are placed on 
system demands, spills and storage levels- 

iii) The flows necessary to meet the levels required hy 
(Eq. 2-3), and at the same time minimise the system's total 
cost of water transport, are determined through the applica- 
tion of an optimisation procedure. 


The mathematical formulation of the directed capaci- 
tated network problem is as follows 


MINIMISE : 


Z 




( 2 . 2 ) 


Subject to 


S q.. . - £ q.. = 0, j = 1, .... n ... (2.3) 

i i 31 

••• 


where 

q. . = Blow from node i to node j for time At( 


C. . == Cost of flow from node i to node j for time 

1 J 

At, 
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Ij. . ~ Loijei' boimd of tb.e a3:*c fj’om node i to 
1 2 

node 3 ; and 

U = U-p-pen iDoxind of the arc from node i to node j 
10 

The first set of constraints (Eq.2.5^ satisfies 
continnlty of mass at all nodes in the network. Eq. 2.4 
describes the unner and lower limits on flow in all arcs in 
the netvTork. The objective fun.ction to be minimised is 
expressed by Eq . 2.2, 


iv) All necessar^r state variables have now been determined, 
and the stat-us of the system at the conclnision of the current 
time step becomes the status at the beginning of the next time 
st ep . 

This procedure is repeated in a step-wise fashion until 
a specified simulation interval has been spanned. 

In SIMIll-II the optimal allocation of netvrork flows is 
accomplished through the application of the " Out-of-Eilter 
Algorithm" (Ford and Fulkerson, 1962). This procedure finds 
the minimum total cost of vrater circulation within the network 


system subject to flow constraints placed on the system arcs. 
If we define the amount of water flowing in arc (ijj), from 
node i to node j, as q. . , and the unit cost of moving water in 
this link as G.., then the algorithm minimizes the objective 

J 

function 1 C q for all i and j in the s^/stem. This is 

X J X J 

accomplished subject to the flow constraints given in 
Eqs. 2.3 and 2.4. The ‘Out-of-Filter Algorithm‘ requires 
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the objective function and all constraints to be linear and 
therefore, SIMYID-II can be considered to be a, linear pro- 
gramming formulation, 

2 . 2.5 Prograrame description 

SIMYLD-II consists of a main programme and eleven 
subroutines, all of which are written in Portran IV progra- 
mming language . Pig. 2. 2 shows the organization of the 
code the subroutine names and calling programme. 

The following is a description of the important 
features of each of the subroutines. 

SIMYID-II (Main Programme) 

The main programme is the control point for call- 
ing subroutines. The Fortran logical unit requirements are 
read in and their values are kept throughout programme execut- 
ion. 

Subroutine ADJUST 

This subroutine is used to adjust the annual de- 
mands in the firm yield calculations. This adjustment is 
based on the greatest shortage incurred during the period 
of operation. When the shortage demand ratio is within the 
user specified tolerance, or the preset value of 10 pe rcent, 
it returns to subroutine OPEATE with input = 0 value. 
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Subroutine CARDS 

Subroutine CARDS reads all input from Cards 
except for the monthly ■variable data( inflo'ws, demands, 
and evaporation rates )• 

Subroutine DATA 1 

This subroutine is called only if variable monthly 
data ( inflows, demands, and evaporation rates) are being 
read from cards. The data cards are read, rearranged, and 
a temporary scratch file is written for use by the programme- 

HTTRY DATA 2 

This entry point within subroutine DATA 1, permits 
the programme to read one year of monthly data during simu- 
lation. The temporary scratch file is written from a previous 
call to DATA 1 or is created in advance. 

MTRY RULE 

This entry point is where the monthly operating 
rule criteria is set- The preselected subsystem of reservoirs 
is analysed to determine if it falls in the average, dry, 
or wet state and the appropriate operating rule is passed 
on to OPRATE. 

Subroutine OPRATE 

This is the major subroutine in the programme 
and is where the yearly and monthly loops are set and all 
calls to" operating parts of the model take place. The arc 
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■bounds and unit flow costs are calculated 3n. this routine. 
This su'broutine also controls the calls to ADJITST and 
DATA 2 for iteration if the model is operating in the firm 
yield mode. Upon return from the network flow algorithm, 
initial arrays and yearly total summaries are set up. Baser 
on an input operation , selected years are passed on to 
OUT 2 for printing. 

S-uhroutine AREA 

This subroutine performs a linear inteipjolation 
to determine reservoir surface area as a function of 
volume. The user is permitted to input up to 18 matched 
points per reservoir to describe the area - capacity re- 
lationship of each reservoir in the system. 

Subroutine OUT 2 

OUT 2 prints detailed monthly system operation fo2 
selected years. At the end of each simulation year, OPMTE 
determines if the year should be printed and issues a call 
to OUT 2. Detailed monthly information for each node and 
link is printed including beginning and ending storages, 
inflows, demand, spills, transfer amounts, etc* 

Subroutine OUT 1 

This subroutine provides a complete printing of 
the input variables that control the simulation. These 
variables have been read previously by subroutine CARDS. 
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S-ulDroutine OUT 3 

Subroutine OUT 3 is an output subroutine called 
from SIMYID-II ( Main Programme ) after the simulation is 
completed. This subroutine prints summaries of the annual 
operation of each node for each year in the simulation 
period. In addition, simulation period totals and maximum 
and average flows in each link are printed. 

Subroutine SETNET 

This subroutine is called to set up the basic net- 
work system. The configuration is determined by the number 
of nodes and. physical links which are joined to mass 
balance nodes by artificial links as described previously. 

Subroutine SUIERKIL 

This subroutine finds the minimum cost flow in th© 
network. As described previously, costs are determined from 
the ranking priorities supplied by the user. The routine 
is called by OPRATE for every month of simulation where 
the bounds and priorities have been assigned. 

Subroutine RIGHT Entiy LEFT 

This subroutine and entry point are intimately 
relate d_^ to the operation of the subroutine SUPERHL. 

2.3 Capabilities and limitations 

SlMYIiD-II is capable of simulating the operation 
of a system of interconnected reservoirs and non-storage 



30 


junctions with a maximum of 30 nodes ( reservoirs and 
non-storage junctions) interconnected hy a maximum of 
45 links. The maximum simulation period, using monthly 
intervals is 30 years. 

Unit identifications have heen purposely deleted 
from output files so that the model can nse acre-feet, 
thousands of acre-feet, metric system units, etc. as units 
of flow and storage as required by the user. 

The model allows for three sets of operating 
policies to guide the model through corresponding hydro- 
logic states. These hydrologic states are deteruined by 
the programme each month from the value of storages and 
inflows in a specified group of reservoirs* 

The model accepts card input for the monthly data 
for a system consisting of upto 20 nodes and for a simula- 
tion period of less than 21 years, for larger systems or 
longer simulation periods, a binary tape must be prepared. 

In some cases, the algorithm, in SIMYLD-II causes 
the simulation to be terminated due to infeasibility in 
solving the network problem. In almost every case, these in- 
feasibilities are caused by in proper specifications by 
the user* The most common are I 

i) Improper system configuration ,* 

ii) The user has not allowed an adequate number of 
spill nodes ,* 
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ill) A minimuin canal capacity is too binding * 

iv) An unregulated inflow occurs where there is no 
possible way to reallocate the water in the system, that is, 
no spill node has been provided, and 

v) Basic data problems 

Most of the above problems, can be readily avoided 
by carefully examining the schematic diagram and verifying , 
that the data are correct* 
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3. IMPIMMTAHON OP SniYID-II PEDGR/MT^IS 

3*1 General 

SIMYIP-II is a fairly lengfhy programme involving 
aronnd 1400 PORTR/GT Cards. Available core memory in IBM 7044 
at I.I.T. Kanpur is less than 32 K. vroi^ds. This is not 
adequate for the original SIMYID-II programme* Two compilers, 
are available for IBM- 7044 viz* (i) WATPOR, (ii) PORTR/iN IV. : 
WATPOR Compiler gives an extensive diagnostic report hut 
memory availability is less and diagnostic checking is slow* 
Due to memory problem, and even for debugging purpose, the 
whole programme could not be run at a time. Main programme 
and each subroutine wei*efed separately on WATPOR compiler 
for correcting the syntax errors. 

3-2 Modification due to Computer System 

SIMYLD— II was developed by Carles D. Puentes of ■, 
the Systems Engineering Division of the Texas Water Develop- 
ment Board on a third generation Computer System. This is 
to be implemented in the IBM 7044-1401 system at I.I.T. Kanpur 
using Portran IV compiler . While correcting the syntax errors ; 
of main and subroutine programmes, problems arose • 

i) 0 ?he notation for continuity card is changed finm 

\ : 

A,B, C etc. to 1,2,3 etc. to suit the computer system I 

ii) Entry statemeni^were used in the original programmes 
e.g. ENTRY DATA 2 and ENTRY EDIE in Subroutine DATA 1 and i 
ENTRY LEPT and ENTRY DUMPO in subroutine RIGHT. Since entry | 
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statements are inadmissible in the IBM 7044"- 1401 system, . 
the following modifications were made to the programme 

a) Since ENTRY BATA 2 and El^TTRY RUIE which occur in 
subroutine DATA 1, are called only one© each and that too in 
subroutine ORRATE, those were incorporated in subroutine 
OPRATE itself 

b) ENTRY BEET and EN TRY DIMPO have been implemented as 
separate subroutines j 

iii) Some read cards had to be changed to srit the computer 
e.g. READCKEN, 11, END = 22) ( TITIE (l), 1=1, 20) was chaged to! 
ASSI® 22 TO IOC 

GAII EYMSET (IOC, IELAG,-7,38) 

READ (E[N,11) ( TITIE (l), I = 1,20) 

iv) After removing the syntax errors, the complete 
programme was fed and it was found to result in large memory 
overflow- The original programme is capable of simulating 
the monthly operation of a system of interconnected reser- 
voirs and non-storage junctionswith a maximum of 30 nodes 
and 45 links over a period of 30 years* In order to reduce 
the memory requirements to the capacity of IBM 7044-1401 
system, it was necessary to reduce the dimension of the 
variables* It was found that a system with 13 nodes and 

20 links can be simulated using monthly intervals over a 
period of 13 years* 
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3-3 Validation of the Original Programme 

The original publication describing SIMYID-II 
(T\4DB 1972) gave the test data for a hypothetical system 
and the programme was run using the test data. The results 
from the simulation run agreed with the results given in the 
original publication and thus the programme implementation 
was validated. 

3*4 Additions and Alterations 

3 . 4.1 Water year 

In the Original programme the water year is counted 
f^uai first January to thirty first December. In this study, 
the water year is modified to begin on June 1st of one year 
and ends cnthe thirty first May of the next year. 

3 . 4.2 Canal system 

SIMYLD-II considers man made and. pumped canals 
which are costly. Hence these are to he avoided and so the 
original programme minimized pumpage cost for canal flow. 
In certain systems e,g.,Bhakra Beas Sy-stem, diversion may 
be preferable beeause of an existing high head reservoir 
and resultant flow through an additional power plant* Further 
more the unit cost of flow in the river reach and canal were 
taken as one in both the cases. In the revised programme, the 
cost of the river reach is kept as one while for the canal 
it is ten and the programme was modified so that it may not 
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minimize canal diversion- This requires identification 

of some river reaches where diversion is "an degi rah le as 
canals and vice versa. 

5 -4-3 System state 

The operation of the system depends upon 
whether storage and inflow in a given month indicate the 
water available to he helow, at or above average conditions 
referred to respectively as dry? average and wet states* fhe 
states are determined as follows * 


A specific group of reservoirs are used to identi- 
fy the system state, 
let 

I +]^ 

= capacity of the i reser*voir 




^i,t 


Y 


i, t+1 


W 


number of reservoirs used in identif^ring 
the state of the system 
the number of reservoirs in the system 
the period under consideration 
end-of-month storage for the i 
reservoir in the t fh time frame 
Unregulated inflows to the i 
reservoir in the ( t +l)fh- time frame 
total storage capacity of reservoirs 
defining the state 



R 


N 

I 

i=l 



11 

+ £ 
i=l 


Y 




36 

( 3 . 1 ) 


let and be fractions of the grbsystem maximum capacity 
used to determine the limits of the hydrologic state with 
X^ _< X^ . Define 

Lg = X^ ¥ and 
= Xg.W 


vrhe'f'e the hydrologic state is determined by 
R < Ig = DRY 

% = AVERiiGS , snd 

R > = ¥E1 

Ju 


Associated with each one of these hydrologic states 
there is a corresponding set of operating rules and ranking 
priorities for meeting demands* 

Based on input parameters supplied by the user for 
the operation of reservoirs and priorities, the programme 
optimises the reservoir releases. 

In SIMYID-II the decision about the states were 
made monthly but in India with a highly seasonal inflow, the 
state of the system does not generally vary in the non- 
monsoon sea^son. Bor example the decision concerning the Bhakra 
Beas System are generally modified month after month ( or 
every ten days ) in the filling season but may be considered 
as constant in the nonfilling season. Hence the programme 
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has heen modified to define and evaluate the state in the 
months of June, July, August and Septenher and use the 
state in September for all subsequent months until next 
June. 

As rule curve values change from month to month, it 
was necessary to define ¥ in tems of desired storage in 
an average year rather than the total storage* Hence 
V/ was redefined as 

F 

¥ = I S. , where S. , is the desired storage 
i=l 

for ith reservoir in t time period for an average year 
as per the rule curve. 

3*4*4 Energy 

SIIIYIE-II prograinme was developed to consider the 
demand fon water* say for irrigation, water supply, etc* 

While irrigation is the major purpose of Beas Sutlej system, 
the energy needs ar© also important and it is necessary to 
keep track of the energy generation from month to month 
as well as seasonal and annual deficits in energy. -Subroutine 
MERGY is added to suit the specific requirement of the 
Beas Sutlej system and using the specific ' Elevation- 
Efficiency - Storage ^ relationship of the system. It may 
he noted that for other systems the subroutine -’■s will 
need modification. * 
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3-5- Details of Resiilts from the Comp-uter Programme 

The o'utp'ut from SIMILD II consists of three 
suhreports as follows * 

Suh-Report 1 

This suh-report contains the information supplied 
hy the user for the simulation and includes the number and 
description of the nodes, links and configuration of the 
system. 

Sub-Report 2 

This sub-report gives details informations the 
nodes and links on the monthly basis and this includes initial 
storage, unregulated inflow, upstream spills, demand, surface 
area, evaporation rates, evporation losses, downstream spills, 
shortages, pumped into, pumped out, system loss, end of 
month storage and operating rules for each storage reservoir, 
demand, storage, and unregulated inflow for each non-storage' 
nodes, and actual flows in each link and the annual average. . 
This is repeated for every year. 

Sub-Report 3 

This sub-report provides a series of summaries 
the period of simulation and also for every year of simulation 
for each node and for each link. 
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5.6 Concl-usion 

It is possilDle to implement the o-^i,ginal SIMYID-II 
in IBM- 7 O 4 A system and adapt the ppoghamme svltahly 
thon^ modifications and additions in onden that 
Bhafcra Beas System can he simulated. Final nestilts 
considered in Chapter-4 indicate that further modi- 
fications should he ro.ade in the programme in order 
that it is more realistic i However the exnerien.ce 
with SIMYlB-II indicates that when a large number of 
computer programmes are available for simulation of 
water resource sxrstems, it may he easier to adapt some 
of the existing programmes than to write a new 
programme particularly when the systems are complex 
and not well understood. 
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4- SIMULATION OP BEAS SUTLEJ SYSTM 

4*1 System Description (Bhalla and Bansals 1975 , ^Mehindiratta. 

and Boon 1.973 a, Hanbans Singla,1964;B'BD05l964 )» 

The Beas Sutlej system (Shown in Pig- 4»l) has "been 
chosen for the study. The river Sutlej, which originates in 
the regions of Mansrover in Tibet, enters Indian territory 
near Shipki and after flowing for a length of about 200 miles, 
it emerges in the plains of Punjab at Bhakra.. Total catchment 
area of Sutlej above Bhakxa Dam is about 21,360 sq^ -miles and of 
this 14,305 sq. miles lie in- Tibet and only 7655 sq-miles lie in 
India. 

The Sutlej catchment is affected by summer 
rainfall as well as winter rainfall. The period of south west 
monsoon rainfall extends from June to September, and winter 
rainfall extends from December to February. Rainfall in the 
catchment varies over the basin with an annual average of 
around 875 mm. Govindsagar, the reservoir formed by the. Bhahra 
dam has a gross capacity of 7.644 m. a. f. and a live storage 
capacity of 5-932 m.a.f. above a dead storage of 1462 ft. It 
covers an area of 41,000 acres. The total runoff at Bhakra 
for a dependable year works out to 11.128 m.a.f. and that for 
a mean year to 13-329 m.a.f. 

Water from Govindsagaji* can be passed through turbines 
of two power houses, one on the right bank and the other on 
the left bark at the foot of the Bhakra dam. Both the power 
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1101:13 63 have 5 tiirhines each. Each of the generators in the 
right bank power house has a maximum capacity of 120 M .lil. where- 
as each of the left bank generators has a maximum capacity of 
90 M.¥. 

About 11 kras. downstream of Bhakra dam is Nangal reser- 
voir formed by the 9 5 ft. high Nangal Dam. Part of the water 
from Nangal is released to Nangal hydel channel with a length 
of 40.07 miles and a carrying capacity of 12,633 cusecs. The 
remainder of the water is released to Sutlej. The Nangal hydel 
channel supplies water to two power houses on its path at 
Ganguwal and Kotla with a total installed capacity of 154 M.W. 
Water from the Nangal hydel channel is then divided, between the 
Bhakra main canal and the Sirhind Canal. 

Downstream of Nangal, there are head works at two 
places on the river Sutlej at Rupar and Harike. At Rupar water 
is diverted to Bist Doab and Sirhind Canals. 

The Be as river takes off from the lower ranges of 
Shiwaliks and joins the river Sutlej at Harike. The total 
length of its course upto its confluence with River Sutlej is 
about 247 miles and the length upto Beas Dam at Pong is 143 
miles. The catchment area of river Beas upto Pong is approxi- 
mately _4, 8 50 sq. miles. The average rainfall in the catchment 
is 1, 778 mm. For a mean year the discharge at Mandi Plain 
varies from 5,328 cusecs minimum in the dry season to 65,350 
cusecs during monsoon, with an annual average runoff of 
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13*01 in.a.f. For a dependable year, the runoff is 10.00 m.a.f. 

Beas project has been undertaken for harnessing the 
water and power resources of the Beas river by means of storage 
and diversion works. It consists of (i) Beas-Sutlej link, 
which comprises a diversion dam at Pandoh across the Beas in 
the Kulu Valley to transfer 3.83 m.a.f. of water to the Bhakra 
reservoir through tunnels and open conduits capable of passing 
a maximum discharge of 7500 cusecs and (ii) Pong dam which 
provides for a storage dam at Pong with a maximum height^of 432 
ft, a gross storage of 6.9 52 m.a.f* and a live storage capacity 
of 5*908 m.a*f. The power plant has 4 units with an installed 
capacity of 60 M.W. each with provision for two additional 
units in future* The water released from Pong dam and utilised 
for generation of power will be used for irrigation through the 
Beas Canal system from the Harike head works. Water from the 
Ravi river is transferred by diversion at Madhopur head works 
through Madhopur Beas link (maximum capacity 10,000 cusecs) to 
the Beas river. This can be diverted at Harike to irrigate 
Beas command. The interconnected system of the Beas 9 Sutlej 
and Ravi rivers is shown in Fig. 4*1* 

4*2 Node link Representation of Bhakra-Beas System 

For the consideration of the study the Beas-Sutlej 
system may be considered to consist of the following ; 
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FIG 41 INTETR CONNECTED SYSTEM OF RlVERSRAVlj BEAS AND SUTLEJ 

(;'BHALLA ANl^ liANSA, ;f75j 
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Two reservoirs at Bhakra (l) and Pong (2) and diversion 
point without storage at Pandoh (3) and five control points 
without storage at Debar (4), Nangal (5), Rupar (6), Harike(7) 
and Mukerian (S). It may be noted that Rangal reservoir has a 
small storage capacity with reference to releases in a month 
and it is used for within the week and within the day balancing 
of demand and inflow. Hence the fluctuation of stora.ge at 
Fangal is. ignored in the study, and Nangal is treated as a 
control point without storage. 

The number within the brackets is the serial number of 
the node with reservoir numbered first and other nodes later. 
The numbering of the river reaches and canals were done initia- 
lly by treating the diversion canal as a canal and river 
reaches as river reaches. However because of the formulation 
of the original programme which minimises canal diversion, this 
formulation resulted in large diversion through Pong and small 
diversion through Debar and Bhakra. 

The inflow in river Beas may be diverted through Beas- 
Sutlej link for augumenting storage at Bhakra and hence to 
generate power at Debar and Bhakra or permit to flow through 
Beas river to Pong, but diversion is always preferable and so 
flow from Pandoh to Pong is desirable only if the inflow at 
Pandoh is larger than the capacity of Beas-Sutlej link canal 
or when there is no storage available in Bhakra. Hence for the 
system model used in the study, Beas-Sutle j link is considered 
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as a river reach with a high priority and Be as river Between 
Pandoh and Pong is considered as a canal section with a low 
priority and is numhered last as shown in Pig- 4*2. 

4-3 Data Used in Simulation Analysis 

Simulation analysis of Bhakra Beas system needs exten- 
sive details of data concerning the reservoir, canals, power 
houses, demands, inflow and priorities. Those have Been obtain- 
ed from Bhakra Beas Management Board (BBMB) and from a number 
of publications (BBDO, 1964 ) and includ.e the water power studi- 
es of BBMB for cycle point 1921-22 to 1959-60. A brief descri- 
ption of the data used in the study are given in the following 
subsection. 

4*3.1 ■ Irrigation demand 

The irrigated area served consists of parts of Punjab, 
Haryana and Rajasthan. The main crops grown in this region 
are bajra, cotton, maize, rice, joT-.ar, sugarcane, oil seeds, 
pulses, potatoes and fodder during rabi season. 

Beas Design Organization (B.D.O.) has estimated the 
demands for various command areas of the Beas-Sutlej system 
for dry, dependable and average years- A part of the demand 
is met from outside the basin including western JairiiUna canal 
and the planned operation of the system for dry, dependable 
and average years are given in Tables 4*1, 4*2 and 4*3. 

The requirements are low in the months from December 
to April except in the latter half of February and in March 
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INFLOW 



FIG. 4.2 node, link REPRESENTATION OF BHAKRa BfeAS SYSTeM 
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when water is required for maturing of rabi crops* The require- 
ments are high in the months from September to NoTember owing 
to water required for maturing of Kharif crops and preparation 
and Sowing of Rabi crops* Rabi crop requirements are less than 
Kharif crop requirements. In the months of May and June, the 
requirements are again high due to preparation and sowing of 
Kharif crops* In the monsoon season also the requirements are 
quite high as the areas where such water is to be utilized, 
have low rainfall. 

In the absence of additional prior information, the 
year is identified in this stufly as dry, dependable or average 
and the net corresponding B .D.O. demands were used in simula- 
tion analysis. 

4*3*2 Power demand 

The releases from the reservoir were originally planned 
to be made mainly in the interest of irrigation. Water is 
released in the interest of power from Bhakra reservoir only. 
These releases were to be in addition to those made .in the 
interest of irrigation. Only 0.154 m.a.f. was earmarked for 
releases to be made to firm up power generation from Bhakra 
for average year while 1.084 m.a.f. and 1.099 m.a.f* was 
planned to be released in the interest of power for dry and 
dependable years respectively. The details of releases in the 
interest of power are also given in Tables 4*1 to 4*3* 
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4*3 *3 « Inflow data 

The Beas-Sntlej system receives inflows from Sutlej 
at Bhakra, from Beas at Bandoh and Pong^and from Ravi-Beas 
diversion and Beas inflows at Mnkerian* Inflow data for river 
Sntlej at Bhakra and for river Beas at Pandoh and Pong are 
available for ten daily perio(3sf«fr a large number of years. 

In this study monthly data for only 13 years from 
June 1961 to May 1974 were used to derive the Sutlej inflow 
at Bhakra, Beas inflow at Pandoh and the additional inflow 
between Pandoh and Pong lumped as Pong inflows. 

Rivers Ravi and Beas are interlinked through M .B .link 
of 10,000 cs capacity. B.D.O. recommended inflow at Mukerian 
through M.B. link for dry and dependable years as 1.442 m.a.f. 
and for the average year as 1.949 m.a.f. Details for monthly 
values are shown in Tables 4*1 to 4.3. 

4*3 •4* Evaporation rates 

Evaporation occurs from the reservoir, channels and 
irrigation areas* Evaporation from irrigated areas is includ- 
ed in irrigation demands, but the losses from reservoir storage 
due to evaporation constitute a consumptive use and are not 
available for other uses. Hence they are to be estimated for 
Bhakra and Pong reservoirs. Mehndiratta and Hoon (1973) 
measured, the pan evaporation at Bhakra from April 1966 to 
March 1971 and estimated the evaporation rates from inflow, 
outflow and storage data. They estimated an average pan 




Maan head (fi) 

2760 - Col 8 - 

fail race f/t'yaifon 


STUDY FOR DRY YE Ah 

















































JUTlCi-JO 







13/77 


733/2 iBB/9 


ttOti ttSZB 


8Z( O SS63 


6630 1 70 SB 


loots fO $! 3 


$383 


4 

IV ^ 


IZ 4^33 12 323 0 i 2 


2 ( OS 4 Zt SI 6 ISIS 


I 3 177 \ 137 01 


/3 

14 



zom 


19643 


7990/ 

64S3 


m 


t3707 


13707 


fSB/S 


131 


9 


73### 


4 

zoo 


400 

> 

ZOO 

, ««•* 

3 00 


$BS 


iOBf3 


7203 


416S 

fss/z 

il3* 

/S36t j 


13034 

ISS4} 

13% SO 

7 375“7 

/S9O0 

/&70S 

Z s 439 

Z7Z04 


- SQO 


*-6 00 


- eoo 


- €00 


~g: 


14-141 set I 


I43S7 


I730S iszet 


/97S2 


z/esf 


Z7S4- 










































































































































+ IS3Z. 10060 3Zil ISSO 
+ IS80 ~ 1 0000 4791 ISH 

■f- to9 ^ 4000 4.S96 I its 

+ /s9 • 3000 orooz /til 


7a 3000 S' 37 IS 13 r 

177 ^OO 4952 isn ^ 

555 " eooo 4-£S6 TsFi 1 

4€3 6000 4/87 I530~i 


m4 6000 3 S 3 7 1566 I Hi 

I >• »* * * > " *" ' "*-, *^’*^ 

660 *6000 tSBl 1536 fl7Z 

577 8000 2 2 9 6 fSOi 1172 

440 tOOOO 1666 7473 /HZ 


6 lOOOO nOO 1662 I 


1736 t4€S-\ niz 


t46Z me 

7662 1179 


3000 


17A. >7 ThO (JB 











































































































































Ht A O 

,..OSS IN ruNNgu 


EPNDABLE 























P.A.W^ 











































PERIOD 




fen o o o o -t .-tt 








































qA .“s 
</% 


T 


si x^* Si 

■< \ i? 

gfS 


iT ■ rv 

■5 M § 
^^i:' -5 



« I .§' 5 j§|| 

SyiAQ-<a»[>S 

iS ie i 


s- ^1 
^7 2 8 


Q> 

\ X I’ 
^ ^.p 
^9 


"f 2 Soq\ 

'^ 2600 i 


I 

■ 5 X 

k> r 

^ ^ h» • 

lA ^ ^ ! 

^ r, ! 

I CN <V ^ > 

“=+^=s«==^ 

30 j 
! 

;?t i5ooo! 


o 
<u 2 

5^ 

y 

\ ^ - 

OS '< 


V 

^ Si 




Kf 

%i .0 
J -K 

>» ^ 
f§^ 


=f” 





__8647!'f^24055] 4 1491 j 15000! 

8 647 1+4597 3 i" 28 50 1 ^ SOOO 
JL3S09jf^553 _ t 4 9! j , 5000! 

!a93l j__64l4 f 128 j '5000; 

! 

■ 1 , f 

■ t 

2 5f 48 Z Al§ 

2 2 1 , 84 ! ,riS^8 9 ___::: !L6li,._J |50.oi_ 

-• *^?4 Z_. 84 8 1050 0. 

:I587!-8995T ~ 5 58' !0500' 


'•■T 
24 26! 

5262 

5748 ^ 


42000! 

; !a93li 

j 6414 

1 



_J;200q’ 
+ fOOO- 

p,.,25f4J^ 
j 22 1841 

-!q903'_ 

-|5589 

i 500^ 

, 1 68 82: 

»-l4t28__ 

4* 5 CD 

f' 

i : 1 5 8 7 

-*899 5 

{ 

}+ 505 

! 

i 

j .34 42 

1 ‘*#404 


!3I2 ' 

J0S71 

1372 

1087 ^ 

1383 

i089 

I3SS 

1092 

!38l 

\ 

( 

I 1094 


4668| i360i 1093 


3 81 Oj^ 
3 242! 


13-^2, 1091 

i330l 1089 


\^ 500| 96! 8 I -610^4 

+ 5001 !70I7 L|373U - 

! - ! i+ 2825 [ + 


40 4 j - _ as 
6!03f4 342 


85 1 
169 


' IDSOo] 

105 00: 

laoGo '!’' 

isoooT 


■■"4' 

. 500l 

500 1 
5 oo ; 


!2 7 i 4 : *-92 80 '! 
19 5021-13672 


, ■^0 2 4 .>■ 25 0 ! 4 
1126 sl^.'^-Sol +• 


0 - 58 .^ 


5? 5 I' 15000 
' 273 ! -SObO! 

5j ' sobo^ 

’ ■■ ' i, 0 - 0 93 1 


3206 i 

1,32 9 1 

1085 

2 8 54 

1322 

1088 

\ 9 89 1 

■ 1303 

10 9 1 

2 1 4 41 

1306 

1085 




155 4 

!. 1 288 ' 

1089 

J2__7 6 

1 12 77 

j . ro92 

9 50 

[ 12.6 0 

! 1087 

5 5 

9 . .50- 


i 1088 


C s ^ ^ .s 


A t ^ Mlih'on '^cre' 


t.a. f. 


§ ij 

M K 



2§.:^: 


^ ^h.S ht; 


: ‘ * i 

3 ?? ^ J 217* 1^0 0*84! 

. ^|_. ^_,L „ .^2 ' ^Al- 0‘ 8 4 1 

I- ^ ^.i..,.. i 28 oT~ 0 * 82 !” 

I 296| 10 ; 


■ ft: ^ 

I 

or 

3vS> 


-4-0 


6j ' 285 i CrS 2 i 375 


} 52 { AtiS j 

ITli' 


2 67| 
251 I 

I 41 T 


244 i 



•234^ 

2 I 2| 

221 i 

— f 


2JI.L.„2a24 __ 0482. 

257j 2671 o-eaf 


Id 24 


249 


2 31 I 236 


: 0*3 4‘ 

■f .._ -i_ 

; o-84i 


sosjS^y^ 

300 /^Q\/ 



i 94I Dec. 



2 38 

2 35 

, 0 84 

57 

22 4 

1 

' 231 i 

j o' 84 i 

{ 59 

202! 

2 13 

10-8 5 

26 I 





2 2ii 2 111 085 i 


189^^65 ■ 0 85_j 
£7 51 . l a z I . 0*8 3 I 


-Jp—-' 


i72 ' l o" "' 


Sane! Acre De--‘f' 


167 


177 *0 81 I 

T6“4"*“’dS*r'" 

■ 

J . 




2 50 yane /' /O 


; m 

'■% i U 


X 

'*^ d. 

^ ,U« 

« «j, v' 


5 UJ V'S 

® S ^ 

h' < 

< X 


<>»mJ liwaM 

1*1 5> 

££ ~ 


cr O 


9 O 94 
2851 


9453 

6255 




2575 


'O’- 154 


M. 4.' F' 


-j J /-.r. z 




4 , 3 COmV D . , 


Q |«-| 


tj" r% 





Qr.SZ“f 

'ut^RoH 


I H i > ^ 

^ O i f- KJ 

2 B ^ 1 S < ££' 

I ^ _ oc h- Cl 

5r X ^^5^3 
p £D P' i X < CC 

o 8 < t 

4 5 


^ r% 

f ^ I 


ICL O 
! a. cc 
5 3 X 
ifi h- 


|mi» 

■i- ^ 

o 9i 

§ iC« 

f Rb 


i& S^jy's I 

s! 


< 

X 




' Q. ^ CL^Li] r. b: 

laRSo^b g 


— ^ ■^<€C 

5 %Q!ia3i 

' “ O a ^>i 

< iSPi 


^ n. 

y. ^ 

I** O 

y %j 
^ >€ 

& I 

O O 

Z o i 

^ 2.1 
3 U 

Q Cl ; 


— Px a. 
&d 5 a, 

*«^klE ^ 

^ "1^ § 
!Gi^<< 

I H Q> ££ U_ 


si '! 2! < 

8 o§ 

^ ' mm* 

Kits 


$ m 5 

O jin o 

O I (V _ i /5 

' H 2 u 

3 20^ 

P lo a 

it 


I if} < 

OlpJ 

3 CC W 

a 

Of Z' D o 

oy a ^ 

s? 9 y 
|g^< 


" 2 ' 5/1 r\ i kl f/'* 

iR? u >/> :» ii 

■ OQ i -3 E Is: U i a 
= : S " 5 : q H 

' 3 H < 


R ' a a 

S ^ b S 

,4y, tn < X 


\j I hi 
, _J i— 


.'JULY '49 981 

i I 

' AUGUST "" 53263 
ITeR ho: 33 944 




SER 2H3Q 1 9246 < 7 500 
OCT. n©878 I 6454 

HOV. ' 6 2CB 307 7 

FeC ' 4? 38 ' 2259 


7500 

SI A S' \ 

420 

i ^ 2 O 

1 325 ; 

1 22493 1 

22493 

23238 ; 

2201 

-f 150 0 

7500 

60 7 63 ,: 

4 20 



1 4 20 

i 325 i 

1 - . ,. . i 

: 224 9 3 j 

22493 . 

23 23 8 

3595 

41500, 

7 500 

4 f 4 44 

540 

i 5 4'^' 

; 3 25 ^ 

' 263 5 5 

1 26355 

2 72 20 

4045 

4 - 500 

7 500 ^ 

32409 i 

i 540 

540 

3 25 ■ 

; 2 63.55 

'. 2 6355 

27220 1 

i 6 1 5 

+ 500 


' ' ' ’ ' i 

1 

i 

( 


I T 

1 

i 



7 500 

” 26 746"' 

i " 54o'“ 

! 540 

: 325' ' 

' 26355 

263 55 

'"'2'7 " 2 ’2 b" ■' 


” ^ 50 0 

64 54 

17 332 ' 

4 80 

4 BO 

; 3 25 

2 66 24 

26624 

2 74 29 j 


1 

3077 

92^5 : 

4 80 

■ 4 80 

j 325 

24043 

24043 j 

‘24'848 ' 


■f'" 

2259 

6997 

, . 336 

, 336 

: 325 

"T6 4'20 

; 6 4 2^>" " 

' « 70-Q 


* + 100 


195 37 
! 8!43 

2-26 75 
2 510 5 

25 203 
2 7429 

24 848 


^ '{ 


JAM. 

FEB. 

MAR. 


432S 
44 CC 
4 786 


2039 

6367 , 

^7C 

1 270 : 

i 325 

2259 

* 66 59 ; 

366 

, 366 j 

-j , , , j 

1 325 

3402 

8188 

4 20 

420 

3 25 

5671 

i 12594 : 

300 

1 300 , 

1 325 


MAY 15472 

— — 

JU-NEHC I 239 iO 
!l-20r3obo7 

2I'3C^ 38045 
M-A.R 


7500 


2 29 72 


7500 • 3 1410 

, i , '. 

7500' -■ 37 507 

-i 

7500 : 45545 


4 80 
480 
48o" 


480 j 325 
480 j 325 
480 j 325 


480 1 325 


' 13261 
I 17708 
: 20056 
: 13088 

I 

"2 606,8 

j 

26499 

26499 


13 261 

^ 1 3 F 56 

200 

■ L3-5S6 9 0 94 

1 7768 

; ■ 1 8399 

; 4* 40'0 

" I- 7 9 9 9 ' 2 8 51 

20056 

20801 

— 4 200 

' 2060! ' 

13088 

13713 

.J . 

— • - 300 

14013 : 5 6 12 


J 


'"'i ' 1 

' t 

t f 

260 68 

V 26873 ~ ^ 

1 — ^ 600. 
r — ■ H . .500 

2 7473' ' — 

... 4 ...... 

26499 

27 304 

' ,: 27904 ; _ 


2 57 5 


. . 

'b*304 


—■I" ' ’•HI 

5 

j 

L— 

f 


i 26499 
■ 7 15' 5*33 


, 26499 
26499 


27304 

27304 


6_p0 

.600’ 


27904 94 53 

"27904 ■ "'6b'’5’6 


O' 503 


OS - ousccs 


•t— — 


M.. A, F 





55 


coefficient of 0.7 for Bhakra reservoir .ifeing meteorological 
data and normal meterological condition, the average monthly 
evaporation rates have Been calculated . The average annual 
evaporation is around 5 ft. (4*80 ft) and the monthly values 
from June to May are 0.73, 0.43, 0.32, 0.39, 0.34, 0.25, 0.17 
0.16, 0.21, 0.37, 0.63, 0.84 ft respectively. These evaporation 
rates were used for Bong reservoir also. 

4 . 3.5 Elevation area capacity curves- 

Elevation area capacity relationships for Bhakra and 
Pong reservoirs were obtained in tabular and graphical form 
from B.B.O. and B.B.M.B. These were defined by sixteen matched 
points per reservoir and are shown in Table 4*4. Eor inter- 
mediate values linear interpolation was adopted* These data 
were used to define minimum storage, the maximum capacity of 
the reservoir and also the rule curve. 

4 . 3.6 Data for energy generation 

Power is generated in the system at Bhakra, Pong, Debar 
and Nangal including G-anguwal and Kotla power houses. RaDSndra 
Eagar and other small systems in ■^r’-ibiilTdidP'S'S' b'/' are ignored in 
the study. 

The energy generated at Wangal is assumed to be cons- 
tant (154 M¥) as in the B.D.O. studies. Debar is a constant 
head power plant and hence the head and efficiency at Debar 
are assumed as constant and energy generation are directly 



TABLE 4.4 ELB7ATI0E-AREA-CATAGITY OE BHAKRA AND POEG 
RESE.WOIRS 
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Bhakra Reservoir Pong Reservoir 



Elevation 

in 

Et. 

Area 

in 

Ac re s 

Capacity 

in 

Acre Pt. 

Elevation 

in 

Pt. 

Area 

in 

Acres 

Capacity 

in 

Acre Pt. 

1 

1400 

9550 

1037156 

126,0 

19182 

1044000 

2 

1420 

9830 

1232568 

1270 

22041 

1250000 

3 

1440 

11620 

1448854 

1280 

25089 

14,84000 

4 

1460 

11880 

1685808 

1290 

28784 

1754000 

5 

1480 

1482.0 

1955734 

1300 

32337 

2060000 

6 

1500 

159 70 

2266902 

1310 

3 5686 

240,2000 

7 

1520 

18820 

2617696 

1320 

38960 

2776000 

8 

1540 

20870 

3017898 

1330 

42251 

3180000 

9 

1560 

23130 

3461558 

1340 

45821 

3620000 

10 

1580 

25370 

3950592 

1350 

49336 

4096000 ' 

11 

1600 

28050 

4489236 

1360 

52792 

4606000 

12 

1620 

30830 

5082934 

1370 

55744 

5150000 

13 

1640 

34000 

5736628 

1380 

58629 

5722000 

14 

1660 

37050 

H53038 

1390 

61612 

6324000 

15 

1680 

40150 

7231460 

1400 

64,404 

6,952000 

16 

1700 

43400 

8074112 

1410 

67048 

7610000 
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related to their discharge in Beas-Sutlej link. 

Bor Bhakra and Pong reservoir, the head on the turbine 
depends upon the reservoir storage and downstream releases and 
the efficiency of the system depends upon the net heads and 
these vary from month to month. B.B.O. has considered the 
details cf planned storage and releases for estimating the 
net head and efficiencies for diy, dependable and average years 
(Tables 4*1 to 4*3) • 

Prom a comparison of the heads, losses and efficiencies 
in diy, dependable and average years in the B.B.O. studies, it 
is concluded that#for the purpose of the study ^head loss may he 
considered as constant equal to 4 ft and 8 ft for Bhakra and 
Pong power plants. Efficiencies may be considered to vary 
piecewise linearly as the function of the net head. [The effi- 
ciency (■'■/) net he ad (h) relationship adopted in this study are 
indicated in Table 4*5. It is possible to use abetter proce- 
dure for estimating the energy generated more accurately. 

4-3-7. link capacities 

The capacity of Beas Sutlej link is assumed to be 7500 
cusecs in BBMB studies* The capacity of the river Sutlej and 
Beas are considered to 50,000 cusecs each. The lower limit 
for link flows is assumed to be zero. 



58 


TABEE 4-5 HEAE-EFFI Cl M Cl KEIAIH ON SHI PS 


Bhakra 

Power Plant 

1 Pong Power Plant 

Head 

Ft • 

Efficiency 

' Head 

Ft. 

Efficiency 

Below 300 

0.80 

Below 170 

0.81 

300 

0.80 

170 

0.81 

325 

0.81 

180 

0.83 

330 

0.82 

190 

0.84 

375 

0.83 

200 

0.85 

400 

0.84 

215 

0.85 

425 

0.85 

241 

0.838 

450 

0.84 

267 

0.826 

475 

0.83 

280 

0.82 

alove 475 

0.83 

above 280 

0.82 
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4*4 Sfudy of Planned Operation 
4- 4-1 Planned operation 

The planned operation for dry, dependable and average 
years are given in Tables 4.1 to 4.5. It is proposed to study 
the planned operation of the system in order that the character- 
istic of the system can be understood and improvement can be 
made in the operation of the system. The following additional 
assumptions were made for studying the planned operation. 

i) Only Bhakra is considered as spill reservoir, and 

ii) The irrigation demands were initially taken as average 
year demands. However because the inflow in large number of 
years were much less than the average demand., if was decided 
to modify the demand in the 5th year as a dry year demand. 

B.B.O. has suggested rule curves for dry, dependable 
and average years while the study required rule curves for dry, 
average and wet states. The rule curve for average and dry states 
were adopted from B-D.O. studies and rule curve for wet states 
was estimated from these two and used in the study. [These are 
indicated in Table 4.6- 

4 . 4 . 2 . Value judgements 

It is necessary to define the relative priorities be- 
tween and among storages and demands. Prom a consideration of 
the characteristics of the system, the following general con- 
clusions were derived [ 
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i) Storage In Bhakra reservoir is much more valuable than 
in Pong reservoir because of its hi^ head and its ability to 
supply irrigation water to all demand nodes,' 

ii) Demand at Nangal is more valuable than that of Rupar 
and Harike because of power' generation benefits,' 

iii) Releases in the interest of power from Bhakra have a 
lower priority' 

iv) Storage is moie valuable in a dry year than in wet or 
average years,' and 

v) Diversion through Debar should have greater priority 
than diversion through Pong. 

Using some priority ranks for average, dry and wet 
years and keeping in mind the above priorities several sets 
of ranks for meeting the irrigation demands were used and, 
tested. Prom the results, the values shown in Table 4-7 were 
adopted for subsequent studies. 

In order to consider the effect of relative priorities 
between the storages and for differemt hydrologic states, 
sixteen sets of priority ranks (Table 4*8) were used along 
with the planned operation and these constitute the simulation 
study of the planned, operation of the system. The value of 
and Xg which define the hydrologic state of the system 
were varied between 0.7 and 1.0 for X^ and 1.0 to 1.2 for X 2 
respectively. The results of simulation include the details 
of monthly, annual and total period summaries of the process. 
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TiffilE 4 • 7 PEIOEIOjr BANKS POR IRRIG-ATIOIT MD POIER 
MAiros 


Hydrolo- Name of Node 

gic State • 

Bhakra Pong Pandoh Debar Nangal Rupar Harike Muke- 

rian 


Average 
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20 

20 
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4.8 PTgQBITY BANKS FOR ROL E GURVIS 
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4.4-3 Results and conclusions 

i) Various values of and X2 to define the hydrologic 
states of the system were tried and finally it was concluded 
that values of X^ = 0*90 and X^ = I-IO was seen reasonahX; for 
the Beas-Sutleg system’ 

ii) Initially Bhakra reservoir alone was taken as spill 
reservoir. Bor some of the cases considered, Pong reservoir 
■became full after about sixth year and Because of the limit 
of the capacity of Beas-Sutlej link and no spill permitted at 
Pong there is no feasible solution and this indicated that 
Pong should also be considered as a spill reservoir* 

iii) The assumed demand for the simulation are much larger 
than the inflows and in large number of years, there is scar- 
city particularly from March to June and these occured at 
Rupar and Harike which had the lowest priority. This indicat- 
ed that in order to keep to deficit within limits, it is 
necessary to reduce the demands in some of the years. The 
inflows in 8th, 10th and 12th years are less hy around five 
million acre ft each than the average year demands^ and so the 
demand in these years are to be reduced to atleast dependable 
year demands,’ 

iv) ‘Variation of the priority ranking between storages 
for different hydrologic states indicated that *. 

a) Priority of storage at Pong should be less than that 


of Bhakra,’ 
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b) Priority for storage at Bhakra should he lower than 
that for demands and can at best be equal to that for the 
lowest priority demand J 

c) Priority for storage in average and wet years should 

be the same. The priority for storages in a dry year may be 
slightly higher in order to avoid the large deficit at the 
end of water year or in the beginning of next water yearj and 
v) It is necessary to modify the rule curves for a better 

operation. 

4.5 Improvements In Reservoir Operation 

The simulation study of planned operation of Beas - 
Sutlej system indicate that the operation of the system can 
be Improved significantly by modifying the rule curves and 
improving the value judgement for demands and storage under 
improved operating rule curve. 

4.5.1 Modification of the rule curve 

The planned operation of B.D.O. adopted in the study 
is indicated in Table 4 * 6 . It shows that planned storage at 
the end of May and June as well as other months are different 
for the three states. The storage at the end of May is in- 
fluenced by the inflows in the previous year and it should 
reflect the desirability of carryover in wet years and deple- 
tion in dry years, but in June, the state depends very much 
on the carry-over of the previous year and the inflow and 



demand in a year may differ very much from the end of the 
year storage. It is hence not possible to adopt different 
initially desirable storages at the end of June. 

Using the general criteria that planned end of June 
storage should be the same in each reservoir for all states’ 
the end of the year storage should reflect the possibilities 
of carry over and depletion^ the reservoir fills up rapidly in 
the months of June to September and depletes from October to 
May or June, several modifications for the desired monthly 
storage levels for the different states in the two reservoirs 
were tried and the following rule curves are considered to be 
an improvement over the planned operation of B.h.O. (Table 4 . 9 )» 
It may be seen that the reservoir fills up quicker in an 
average year than in dry year and in turn it fills up qui.ckei* 
in a wet year than in an average year. The end of year storage 
in a dry year is less than that in an average year which in 
turn is less than that in a wet year. The general trend of 
the rule curve is similar to that of the planned operation 
except for some further modification in the end of month 
storages for August and April. 

4.5*2 Ranking of priorities 

Simulation of the planned studies as well as the 
improved operation indicate diversion through Pong from Harike 
while storage and capacity were available through Debar and 
Bhakra. In order to eliminate this problem, the priority for 
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demand at Harike was reduced by changing the rank from 16 
to 25- 

In order to study the priority in meeting the demand at 
Bhakra exclnsitelyrr in the interest of power^ four simulation 
runs were done with changing demands as indicated in Table 
4*10 and priorities as indicated in Table 4*11* 

4 , 5.3 Evaluation of system performance 

The planned operation as per B.I.O. meets the specified 
irrigation and power demands. The planned energ7 generation in 
the components as well as the total system are taken from 
Tables 4*1 to 4*3 and are shown in Table 4»12. 

Energy generated as well as the irrigation deficits 
vary over the months and over the years. Eor the purpose of 
the study the following details of the perfoimance of the 
system are considered important (Tables 4*13 to 4*16) . 

i) Total annual and monthly irrigation deficits^ 

ii) Power generation at (.Bhakra and Debar), Pong, and the 
total power generation of the .system in Kharif season and non- 
Idiarlf season and water year respectively J 

iii) The surplus or deficit of power in each of the sub- 
systems in Kharif season, non-kharif season and water year 
respectively, in comparison to planned values* 

iv) An earlier study (Rao, 1976) of the system indicated 
that • 
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(Values as specified by B.D.O. 



TABLE 4.11 MOLIEIED PRIOKETY MKS 
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1) ConjunctiTe use of surface and ground water resources 
is necessary,* 

2) Since excess energy is awailalle in Kharif season, 
this may le used for pumping ground water to meet part of the 
irrigation demand and to this extent, the surface water 
releases may he reduced. The earlier study indicated that in 
Bhakra reservoir in Eharif season about 1500 acre ft. per 
month of water is required to generate 1 MW of power and this 
will pump 6000 acre ft. of ground water. Hence 1 MW of excess 
enerigy indicates a potential for saving 4500 acre ft of water 
from Bhakra releases and thus augmenting storage J 

3) The amount of water saved indicated in column 25 of 
Tables 4*13 to 4*16 increases the head and energy generated 
in the subsequent period and this is not taken into account 
in this study. 

Yet this amount of water is available for irrigation 
and this may meet partly or completely the irrigation deficits 
in the water year. Net water deficit jif any jin the water year 
is indicated in column 26 of Tables 4*15 to 4*16. 

4) This deficit may be met by diverting energy from out- 
side the system and the amount of energy in MW months required 
for meeting these deficits are indicated in column 27 of 
Tables 4.15 to 4 •16. 
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4* 5*4 Discussion of results 

A compafisan of the generated energy in simulation 
runs with planned demand indicates that surplus energy is 
generated even in a dry year during the Kharif period and 
there is adoficiij nfenergy during rabi season in average and 
dry years. Rao (1976) has suggested the use of surplus 
energy in Kharif to pump ground water by tube wells and hence 
reduce- the surface water releases in Kharif for meeting 
irrigation demands. This will result in a higher storage 
during and at the end of filling period. It will also lead 
the availability of larger amount of water and energy in the 
depletion period. 

Results of simulation run I indicate that irrigation 
and energy deficits are eliminated in a number of years and 
are otherwise reduced in other years. For example the irriga- 
tion deficits in years 2nd, 6th, Idtth- and 13th are eliminated 

h , 

and the irrigation deficits in years 3, 7, 9 and ip are 
reduced. It is also seen that RO to -24 months of . energy 
will be required to eliminate the deficits in years 7 and 10 
respectively. The ener.ey requirements in years 31*^ and 9th 
are of a much higher order, viz, 136 and 483 MW months. If 
energy of these magnitudes is available from outside the 
system, it is possible to eliminate irrigation deficits even 
in these years provided that adequate tube wells are also 
available. 
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Sim'ula'tion Eun II assnines a higher priority for main- 
taining storage in dpy years in comparision to vet years. The 
results indicate that the operation is similar upto 6th year 
and then differs from 7th year onwards, while the rariations 
are erratio from 8th to 11th year. The energy deficit is 
much larger in the 12th year in Kharif period and the energy 
surplus is smaller in the 13th year in simulation Bun II in 
comparison to simulation Bun I. Furthermore the amount of 
energy to he imported from outside for meeting the irrigation 
deficits 'igo also larger in 7th, 9th, 10th and 12th years in 
simulation Run II. It is also seen that surpluses and defi- 
cits of energy are generally smaller in simulation Run II 
than those in simulation Run I. 

Simulation Run III eliminates release from Bhakra on 
consideration of power. Simulation Run IV considers power 
demand at Bhakra only for the dry year 5th with the lowest 
priority. The results indicate that - 

i) The fluctuation in energy is larger for Run III than 
Run IV; 

ii) There is energy deficit in the 5th year in rahi 
season in Run III even though there is a surplus in an entire 
year; 

iii) The energy generation in Run III is smaller than that 
of Run IV hy 126 MW months,* 

iv) There is a large system loss in month 3rd of 6th 
year for Run III which is ah sent :r:i.n3 Run iv; and 
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v) The irrigation deficit for 6th year is smaller for 

Run III than for Run IV and comparable in other years. 

From the above results it is generally seen that 
energy releases from Bhakra need not be specified. Depending 
upon the relative value of irrigation deficit versus energy, 
at .best those may be specified for dry years. In the absence 
of such value systems, it is inferred that Simulation Run IV 
is preferable . 

4*6 Conclusions 

The study indicates the following t 

i) Rule curves and priorities aBo given in Tables 4*9 
and 4 •11,* 

ii) Irrigation demands and priorities a'sr- Indicated in 
Tables 4*1 to 4-3 and 4*7,' and 

iii) Power demands and priorities as indicated for Simula- 
tion Run IV in Tables 4*10 and 4.11 lead to an improved opera- 
tion of the system. 

Further improvement of the system performance is 
possible by modification of the simulation programme in the 
following respects '. 

i) Incorporation of better algorithm for estimating the 
energy generation in the system* 

ii) In some time periods diversion through- Tong are- indicat' 
ed while capacity through Debar and at Bhakra are available. 
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In these cases it is necessary to divert flow throng Debar, 
even though the out of Kilter alfaa'ithin indicates diversion 
throu^ i’ong,* 

iii) Incorporation of subroutine to determine energy sur- 
pluses if any during kharif period, utilization of available 
tube wells for pumping ground water during this period to the 
extent necessary in conjunction with surface water releases, 
and modifications in SIMYID II programme for correcting^ 
releases and storages appropriately; and 

iv) Incorporation of flow and demand forecasting as well as 
suitable flexible control rules for the rule curves. 



® results for SIMULATIO t^ rUI, fjo. T 

-^' 1 ” nriTif, mWwSfertwv ' 

'■ ^ POWER *GEKERATED 


Vj- [TOTAL ^ IBRlGi 
[ S' .tNFLOW l^ui DEI* 

Tt* 1 |H 

Iff' hm 

I* [ ;X Pt, 


TOTAL p 5 tFi(ltc;ATION IlIRtGATION DEFlCtT 
fHFLOW iQ tii DEMAMO July TO JUNE 


1 I \^O 770 :^m^Bt 


■ 4 . ' 

--t- 


AMOUNT m I PETAIUJ 
ACHE FT j 

.....^ , j 


^ POWER 'gEWEMTED -—.^ Mg SURPLUSfr) OR DEFICIT C-) m 

^ ^ ^ f 

POWER IN POWER IN TOTAL SURPLUSOR DEFICIT SURPLUSOR DEFICIT SimilSOR DEFICIT acpe ft water [iiwioationdeficit 'ADc:Tiri^r~*^ 

KHARIF non KH/^IF POWER IN KHARIF IN NON KHARIF INWHOLEYEAR that c*n se aftej. usins POWERliEOumED 

^H^PtA PONG Total ^<^<yt|pONG TOTAtr BHAFfTA;] PONG TOTAL BHAKPAjpONG TOTAL SHAKRAPONG TOTAL BHiMCHA pONG TOTAL STORED SY$Of»»LUS SURPLUS ENERGY TO MEET DEFICIT 
&EHAR I ^±.. • •+• + + ENERGY ^ ^ 


POWER IN 
KHARIF 


0a^AR 


m^AR 


4r 

DEHAR 


total bhakra pong 

TOTAL 

SHAKRAPONG 

•+* 


+ 

DEHAR 


DEHAR 

15 16 17 “ 


19 20 

/7315 +4 64 - 19 

444$ 

455 +354 


+ 

DEHAR 


ENERGY 
IN KHARIF 


r ACRB'^ 
\ FT. / 


R'i—— [-t — LJ — : 1 I l_ i_L.Lh 

«„»«• AHv F,i „„ „.,L„ L 7 . 7 : 

'■■'■ '■• ' llrl ‘ , ' ■ i . 1 


M/i^ 


!■ 

1-^4 

|r^ 


'1 — "i ^ — 

Z 7 ^ 694 « \| 6 IS 4 ^^ i 258 Uw U 2 ...> 6 > L 7 - ^ 9 - 

- I 'r ■ ■ ■'-■1 A. I , .1 7 ^ . . 


j.. 

f't 4 f - 37 ^ Sofooo 



0/4^4- 7>t. 


Nl^ 6^77 UZ 7 


^424 hf) 7 ; 


lZ 8 ipn 2444 yyoBS 


4 756 


* / ^ 75775 ^ I h 


jj?f,a:64447.f 5, .j | JOiZ& fa^/^ 

( -i , ■■-":■ '4’.' ■ ■ , 


y^7 j • 44S4 : .;?07 '1 5414 |4r64-[|«T4 f n786 


S6<7 


+- 24 ZK 47 / I+ 5 Z 3 (-+ 7 ^ 41 . 74-840 00 


^ 3 \'^ 7 f 7 oo& 


I 1 tiSf>H4^Z I ^7J7Si ■ ■ ;] C/w| Siyf 

t:- . ■ ■* I " h ■ ■ r ■■■• ■ . /- 


j^ 7 S 5 oa 76 , /v//. 


.-N/L ■■ 1 4^-? 3? 


j 7 >© f. 1 | 6 & 92 i 477 i -h I 4 

••: : j I P'“ i 

76 S 5 I >^ 7 | Moo \^ss 42 . 

' ' 'I 


7 ;74 0004^ ' Wf T,: iZyS^ | 3* 7© 4^ 5" ? L ^ ' 11x40 

: j I ■ , I irj/SA.sr/c^j^y) ' 

u |__„,_j — 1..:. ' ; ' . . ■ ; . 

-) ' 2 )/?’jO' • ^f 9 S'o^ 7 B j S^ 77 ^'‘ 4 - ■'■■ ” '•■” 

ill ! I I 


5274| S5/8 |4:?3 170/9 h/^67 f578 /5793Uz4.y -^6 


f [ ^rT *~1 "" 1 , f I } I '1 

^ 7 oo> 2 ^|'^ j^ 7 S 6 $m 2 ’'jf 2 ,ir 9 f^ * 14.^0 $' 27 o [5951 8SG \ 7 e 7 s I2Ml 46 / 4 ^ — 77 j-F )99 4 ^ 7 ^ 

^f^^3>CTumi j 1 I - I I L 


494 I AO ^000 I {ob-rSt 


T 9 ? 9 |- 43 S j 43 o 4 j 


I 

0 / 4 SOO [ 124 ^Cf 4 .‘y^y 

* 

3-/ r-^So ! 4- <53 900 

-■V. ...,P|-- 

" — ^• ' — r~ ■'■■■ — •■' 




Nll^ 59201 S 3 o 7/70 


94" I 

'■S/;- '■ '■*'1 . 

— 


^717-74 


H362 


? 5568 'S 


WH7 ; 778^94^ /5o47Xif N^y 677/ Uozs 6?566 j S52Z 969 75^9 »225? /954 I607S 4 9 6 1 -f 4 


27 

N/i- 

/ 


Nii- 


/V? t 














































■ 




TA»i.e i . 16 jesuinv.,, fM ..simulation ._RyN ^,Njj,jS: 

< —power" generated — 

irrigation irrigation deficit power in power in 

DEM<\NO JULY TO JUNE I^HARIF NON KHARIF 


SULTt, FOR S 


,.v 


I'tOTAL la liRRKSATtC 

(INFLOW 1 3 ! DEMAND 

I O 

; T, 1 S §- ' FT.'-- . 


KHARir NON KHARIF 


TOTAL 

POWER 


— -—-SURPLUS 05 OR DEFIC — 

SU^F^S ORt^rCIT&Pm^^^^ t|aCI?E ft WMEP ■ TlRRI^^nON deficit jADOmONAL 


|^©UNT m oetails 
; u^m ■ FT. 


IN KHAmr IN NON KHARIF IN WHOLE YEAR power requircp j 

■ • ■■. . ,-, 1 -■■ ■'■A .- ^^ STOPEDSYSURPLUS , SURPLUS ENERGY }tO MEET DEFICfT • 


0 r)#f» 


^ ... ■ ^1 __ I ■ - I — .r ,.,.1 I- r- . ■ --.r -7 -- I ..JI - i -m— -r'-ri - .i - i,-.'.i':; ii ■ r^- -- 1 rmn -1 i .-, i - i n , — rr—p i' m- ii r ti. -Tr-'-f-nr-tT www i miim ii . 

PON® TOm iHAKJM;f»0M6 . TOTAL TOTAL JW<Al{«A i PON® TOTAL|(!>HA»<t«A PONG TOTAL A PONG TOTAL ENERGY 

Fvff'iJAQ ' . 0^"^ HAR^ i — .iF . ... D£l4Afi 


TOTAL 

iHAK,AA-J 

DEHAR 

PO M® . 

TOTAL 

|BHAIL'RA| 

j.'&SwR 1 

k 

PONG 

9 

\d ■ H 

ram 

2..., j 

■13 ' 

t4 


4- 

(DEMAR 


4 

Qlim 


m KHARIF 


WET .a 7 $. 494 «t Ht-L 


UTjaa 


16 "j 

1 17'. 

i 

! 

. ! 

1 

: - i 9 

[ 

1 

i 

1 


23 124 


5 -I a 


S2 +315 +767 2002SOO 

i 


.gj”.: 'bs* 0 «« »S a 7 t *»4 5^2 . ■'r 73 ;t 7 t 


,:': 6,#3 2 , laT'i, bs?# [aosp .’i'-s'T-i. 


4 so 


2 St 1-411 


.231 +. 169 2 <; 1 0000 




'3 ■ b6f 38606 M^j0>4si I T6T 54 2 ,4^ 


I ■ ;•. t 


3 1238 16 271 f6t4t | 768 {7994 H 24 M 12006 M626SI+I6I i + 17 S+17 


NIL 


4 ; iawsis-tt WET-'; 27 S$ 74 fl 

'•,. ■ !. ' . ' 

■'''' ■ : ' ' ■ t ^ 

[ DRY ’!.t 67 5 -r-|S 6 

I ^ I ' ^■'' ’ '/ 

^]^L',r' .i:"' M II ' H 'I I P W w ...... ll■, > 4^w 

I $ |a47iS38.f-, WET:l2,7'5A9 4s'2, \ 2f'«9 62'4 


f lit 7 f 82 74 1630 2.: If,* 2 1 8 772 ;l'2759.j 2439 1 1704414- 27 5 J —94 ^ +»' 8 L ; 4-2 I4 1 •+« 5 3 |+-3 67 1+4 89 !•+ S ? 1+S46 

•■ r' . 1 - ■ • ' -i " I ! ' 

^ ; j ■ . . - ,r„_, i. _ . ; i, 


6928/14801. f 34( 1S92C.99I6 ; JC84 |t2F48 4-443 : + t70 4-613 ; 4 - SC , 4-163' j 4-243 | 4-523 


i * * 

4 

i. • 

’ '>■ I 


[758 5 00 


88115233 i 889 [7200 ll 233U • '. 4-768 - -93 5 j -2 70 -1 205 j +- 6t 


417 13 54 600 O 


7 127359591 fWET * 2 2. ' 5117 5 1 


617 9 '1234 I PI 83 164 64 5 8 67 84© 9 I2S43 : 2101 1165 92 -+77 -413 [-4 90 14.296 j— 292 I -4 4 f .f373 


I 8 ;2J 858098 DETIJ 2 <550876 


iS823 !879 j 7472 59744 [22f TC 43 1 H 567 IIOO 14515 -4-3 


4-3J1. '-i-Sa? 1 “43 7 i+.(5a :-f89j? 


124800488 [WET 1278694 S t 8684721^ 


7028 'II6TO • 1589 ! S3C7 1 4-240 1 — 6 2 i+tP6 


- 1 ' 37 T:- 


lo hit^lOZS 'D'FP. } 7 "ISS 0876 


\ s ’ ; 

‘ NMtl 


o* 1^1 7rf 4-120 


I 6*5 5 0 '.o 


003020 WET [2 7 8 69452 ■ 779 9,7 


« 12387 


+>363 I 


I 7 23 50 O 


f;«.; 1 4% 


2 [23 40657®^} pEp, j ? I ® 50 87 6 


7120 15775 ! 333 7 186 11295 11 6'3 


1 5047 3 9 


7509 12331 i 1934 


2289500 


f I POWER 


;(acr£ 

\ 


6 


NIL 


i +J 41 |_074 -‘233 80 IOCO 1 814422 


106751 


i^fo-f-ccc 1 NIL 


NIL 


NIL 




’ s‘>.. 
























































































































81 


5. SUmARY Ml) CON CUT SIGNS 
5-1' Summary 

Water resources systems are generally large and 
complex- They consist of multiple units and serve multiple 
purposes- The multiple purposes are not wholly complementary- 
Several general or prohlem specific simulation models and 
computer programmes have been developed for analysing water 
resources systems SIMYID II is a computer programme 
developed hy Texas Water Development Board for simulating 
the hydrologic operation of a system of interconnected 
reservoirs within a hasin or a multihasin water resources 
system. The study consists in implementing the SIMYID II 
programme^ validating it with available data and adapting 
it for the operation of Bhakra Beas system. 

The original SIMYLD II programme was implemented 
in IBM 7044-1401 system at I.I.T. Kanpur and was adapted 
to meet the requirements of Bhakra Beas system. 

Using 13 years of historical data, the Bhakra 
Beas sfystem was simulated using the modified SIMYLD II 
model- The criteria for defining wet, average and dry years, 
the rule curves for operation in wet .average and dry years 
and the relative weightages for meeting different demands 
and for maintaining the rule curves are derived from simu- 
lation analysis- Results indicate that by using these 
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criteria the henefit from the operation of the system can 
he greatly increased, further improvement of the model is 
also possible. 

5-2. Conclnsions 

The following conclnsions emerge from the stndy *. 

i) The original SIMYIh II programme can he implemented 
in IBM 7044—1401 system with minor modifications dne to 
computer system and hy limiting the dimensions of the problem 
to consider 13 nodes, 20 links and 13 years of monthly data- 

ii) The programme needed additions and alterations 
a) "because of difference in the water year h) because of 
gravity canal rather than pumped canal c) because of differ- 
ing climate and the corresponding need for defining the 
state of the system to suit the climate and d) incorporation 
of energy developed from the system. 

iii) Simulation of Bhakra Beas system indicates that 
it is possible to derive improved rule curves for the 
operation of the system in wet, average and dry years and 
the corresponding priorities between demands and maintaining 
the rule curves. 

iv) Differing rule curves for wet, average and dry 
years incorporate the utility of carry over from wet years 
to dry years. 
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■v) ’’'i/hen a large number of computer programmes are 

available for simulation of water resources systems, it 
may be easier to adb-'pt an existing programme than to write 
a new programme . 

vi) ’ further modifications for optimization of the 
system and short term planning are possible. 

5.3 Suggestions for Future Study 

i) The present study is limited to the consideration 
of simulation of 13 years operation of Bhakra Beas system. 
Study of a longer period incorporating heuristic criteria 
to' avoid diversion through Pong when diversion through 
Debar is possible, and conjunctive integrated operation of 
the surface water-ground water system and with the addi- 
tional details, will lead to better operation of the system. 

ii) Simulation analysis using generated data for 
stream flow and irrigation and energy requirements, as part 
of a more comprehensive system planning study is also 


needed. 
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